Chapter 36 



Chaotic multiscattering 



(A. Wirzba and P. Cvitanovic) 

WE DISCUSS HERE the semiclassics of scattering in open systems with a finite 
number of non-overiapping finite scattering regions. Why is this inter- 
esting at all? The semiclassics of scattering systems has five advantages 
compared to the bound-state problems such as the helium quantization discussed 
in chapter 37. 

• For bound-state problem the semiclassical approximation does not respect 
quantum-mechanical unitarity, and the semi-classical eigenenergies are not 
real. Here we construct a manifestly unitary semiclassical scattering matrix. 

• The Weyl-term contributions decouple from the multi-scattering system. 

• The close relation to the classical escape processes discussed in chapter 1 . 

• For scattering systems the derivation of cycle expansions is more direct and 
controlled than in the bound-state case: the semiclassical cycle expansion 
is the saddle point approximation to the cumulant expansion of the determi- 
nant of the exact quantum-mechanical multi-scattering matrix. 

• The region of convergence of the semiclassical spectral function is larger 
than is the case for the bound-state case. 

We start by a brief review of the elastic scattering of a point particle from finite 
collection of non-overlapping scattering regions in terms of the standard textbook 
scattering theory, and then develop the semiclassical scattering trace formulas and 
spectral determinants for scattering off N disks in a plane. 



685 



CHAPTER 36. CHAOTIC MULTISCATTERING 



686 



36.1 Quantum mechanical scattering matrix 

We now specialize to the elastic scattering of a point particle from finite collection 
of A'^ non-overlapping reflecting disks in a 2-dimensional plane. As the point par- 
ticle moves freely between the static scatterers, the time-independent Schrodinger 
equation outside the scattering regions is the Helmholtz equation: 

{^^l + 1?^ tf/{f) = , r outside the scattering regions. (36.1) 

Here ifr{r) is the wave function of the point particle at spatial position f and E = 
fp-p- jlm is its energy written in terms of its mass m and the wave vector ^of the 
incident wave. For reflecting wall billiards the scattering problem is a boundary 
value problem with Dirichlet boundary conditions: 

ij/if) - 0, r on the billiard perimeter (36.2) 

As usual for scattering problems, we expand the wave function tp{f) in the 
(2-dimensional) angular momentum eigenfunctions basis 



m = £ ^iine-"'"' , (36.3) 

in=—oo 

where k and O^t are the length and angle of the wave vector, respectively. A plane 
wave in two dimensions expaned in the angular momentum basis is 

CXD 

/■f ^ ^ikrM-<i>,) ^ ^ J^{kr)e""^'^^-'^'^ , (36.4) 



where r and denote the distance and angle of the spatial vector r as measured 
in the global 2-dimensional coordinate system. 

The mth angular component J,nikr)e""^' of a plane wave is split into a super- 
position of incoming and outgoing 2-dimensional spherical waves by decompos- 
ing the ordinary Bessel function Jmiz) into the sum 

Jm(z) = ^ {hI^\z) + hI^\z)) (36.5) 



of the Hankel functions hII\z) and H^miz) of the first and second kind. For |z| » 1 
the Hankel functions behave asymptotically as: 



fi\l\z) ~ /Ag-''fe-f'n-|) incoming, 
V nz 

Hil\z) ~ /— e+'('-f'"-f) outgoing. (36.6) 
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Thus for r — > oo and k fixed, the with angular component J,„{kr)e ' of the plane 
wave can be written as superposition of incoming and outgoing 2-dimensional 
spherical waves: 



^|2nkr 



(36.7) 



In terms of the asymptotic (angular- momentum) components i/'^, of the wave 
function tp{f), the scattering matrix (35.3) is defined as 

^ oo 

~ y + S^^./'^^-'^-'-m . (36.8) 

The matrix element Smm' describes the scattering of an incoming wave with an- 
gular momentum m into an outgoing wave with angular momentum m' . If there 
are no scatterers, then S = 1 and the asymptotic expression of the plane wave e'^''' 
in two dimensions is recovered from ip{f). 



36.1.1 1-disk scattering matrix 

In general, S is nondiagonal and nonsepai^able. An exception is the 1-disk scat- 
tered If the origin of the coordinate system is placed at the center of the disk, by 
(36.5) the with angular component of the time-independent scattering wave func- 
tion is a superposition of incoming and outgoing 2-dimensional spherical waves 

exercise 35.2 

- [j,„{kr) - ^T,„,nHl^\kr)y'""^^ . 
The vanishing (36.2) of the wave function on the disk perimeter 



= Jm(ka) - -T,nmH)n\ka) 
yields the 1-disk scattering matrix in analytic form: 
( 2J^{ka,) \ Hf\ka,) 

Km'^k) = . ^'nm' = - , (36.9) 

V H]n'{ka,)) H;„'{ka,) 

where a = is radius of the disk and the suffix ' indicates that we are dealing 
with a disk whose label is s. We shall derive a semiclassical approximation to this 
1-disk S-matrix in sect. 36.3. 
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36.1.2 Multi-scattering matrix 

Consider next a scattering region consisting of A'^ non-overlapping disks labeled 
s e {1,2, - ■ ■ ,N], following the notational conventions of sect. 11.6. The strategy 
is to construct the full T-matrix (35.3) from the exact 1-disk scattering matrix 
(36.9) by a succession of coordinate rotations and translations such that at each 
step the coordinate system is centered at the origin of a disk. Then the T-matrix 
in Snim' = Smm' " ' T^m' Can be spUt into a product over three kinds of matrices, 



T,„,Ak)-2 Z cJ/^)M-'(^);;>;;„/fc). 



=1 



The outgoing spherical wave scattered by the disk s is obtained by shifting the 
global coordinates origin distance Rs to the center of the disk s, and measuring 
the angle with respect to direction k of the outgoing spherical wave. As in 
(36.9), the matrix C* takes form 



C ; = ^ ^'-~f^^\ i'n^. . (36.10) 



If we now describe the ingoing spherical wave in the disk s' coordinate frame by 
the matrix 

D;^;,,^, = -na,>Jn,-i, (kRs'Vi, (ka,>)e-''"''^^' , (36. 11) 

and apply the Bessel function addition theorem 



Jmiy + z) ^ ^ Jm-e{y)J[{z), 



[=—oo 



we recover the T-matrix (36.9) for the single disk s = s' , M = I scattering. The 
Bessel function sum is a statement of the completness of the spherical wave basis; 
as we shift the origin from the disk s to the disk s' by distance Rs> , we have to 
reexpand all basis functions in the new coordinate frame. 

The labels m and m' refer to the angular momentum quantum numbers of the 
ingoing and outgoing waves in the global coordinate system, and , Is' refer to the 
(angular momentum) basis fixed at the sth and ^'th scatterer, respectively. Thus, 

and depend on the origin and orientation of the global coordinate system 
of the 2-dimensional plane as well as on the internal coordinates of the scatterers. 
As they can be made separable in the scatterer label s, they describe the single 
scatterer aspects of what, in general, is a multi-scattering problem. 
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Figure 36.1: Global and local coordinates for a gen- 
eral 3-disk problem. 



The matrix M is called the multi-scattering matrix. If the scattering problem 
consists only of one scatterer, M is simply the unit matrix M^] = cJ" 5/ / ,. 
For scattering from more than one scatterer we separate out a "single traversal" 
matrix A which transports the scattered wave from a scattering region Ms to the 
scattering region Ms', 

Ki=^''^'^^^'-<i.- (36.12) 



The matrix A**' reads: 

Aff - -(1 - 6''') — '^''^"'^ H^^\ . (36.13) 

as' Hy\kas') ' 

Here, a., is the radius of the i'th disk. Rs and are the distance and angle, 
respectively, of the ray from the origin in the 2-dimensional plane to the center of 
disk s as measured in the global coordinate system. Furthermore, Rss' - Rs's is 
the sepai^ation between the centers of the sth and s'th disk and Os's of the ray from 
the center of disk s to the center of disk s' as measured in the local (body-fixed) 
coordinate system of disk s (see figure 36.1). 

Expanded as a geometrical series about the unit matrix 1, the inverse matrix 
generates a multi-scattering series in powers of the single-traversal matrix A. 
All genuine multi-scattering dynamics is contained in the matrix A; by construc- 
tion A vanishes for a single-scatterer system. 



36.2 A/^-scatterer spectral determinant 



In the following we limit ourselves to a study of the spectral properties of the S- 
matrix: resonances, time delays and phase shifts. The resonances are given by the 
poles of the S-matrix in the lower complex wave number {k) plane; more precisely, 
by the poles of the S on the second Riemann sheet of the complex energy plane. 
As the S-matrix is unitary, it is also natural to focus on its total phase shift r]{k) 
defined by detS - exp^"'^'"'\ The time-delay is proportional to the derivative of 
the phase shift with respect to the wave number k. 
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As we are only interested in spectral properties of the scattering problem, it 
suffices to study det S. This determinant is basis and coordinate-system indepen- 
dent, whereas the S-matrix itself depends on the global coordinate system and on 
the choice of basis for the point particle wave function. 

As the S-matrix is, in general, an infinite dimensional matrix, it is not clear 
whether the corresponding determinant exists at all. If T-matiix is trace-class, the 
determinant does exist. What does this mean? 



36.2.1 Trace-class operators 

An operator (an infinite-dimensional matrix) is called trace-class if and only if, 
for any choice of orthonormal basis, the sum of the diagonal matrix elements 
converges absolutely; it is called "Hilbert-Schmidt," if the sum of the absolute 
squared diagonal matrix elements converges. Once an operator is diagnosed as 
trace-class, we are allowed to manipulate it as we manipulate finite-dimensional 
matrices. We review the theory of trace-class operators in appendix K; here we 
will assume that the T-matrix (35.3) is trace-class, and draw the conlusions. 

If A is trace-class, the determinant det (1 - zA), as defined by the cumulant 
expansion, exists and is an entire function of z. Furthermore, the determinant is 
invariant under any unitary transformation. 

The cumulant expansion is the analytical continuation (as Taylor expansion in 
the book-keeping variable z) of the determinant 



det (1 - zA) = exp[tr ln(l - zA)] = exp 



V n = l 



That means 



det(l-zA) := J]z'"e«,(A), 



m=0 



(36.14) 



where the cumulants Qmi^) satisfy the Plemelj -Smithies recursion formula (K. 19), 
a generalization of Newton's formula to determinants of infinite-dimensional ma- 
trices. 



eo(A) 

QmiA) 



^ m 



> 0,„-y(A)tr(AO form>l. 



(36.15) 



7=1 



in terms of cumulants of order n < m and traces of order n < m. Because of the 
trace-class property of A, all cumulants and traces exist separately. 
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For the general case of < c» non-overlapping scatterers, the T-matrix can be 
shown to be trace-class, so the determinant of the S-matrix is well defined. What 
does trace-class property mean for the corresponding matrices C*, D'* and A*'* ? 
Manipulating the operators as though they were finite matrices, we can perform 
the following transformations: 



detS - det (l - /CM^^D) 

= Det (l - /M^'DC) = Det (M'^CM - /DC)) 
Det(M-/DC) 
Det(M) 



(36.16) 



In the first line of (36.16) the determinant is taken over small ( (the angular mo- 
mentum with respect to the global system). In the remainder of (36.16) the deter- 
minant is evaluated over the multiple indices = (s, 1^). In order to signal this 
difference we use the following notation: det . . . and tr . . . refer to the \£) space, 
Det . . . and Tr . . . refer to the multiple index space. The matrices in the multiple 
index space are expanded in the complete basis (IL J) = {Is,^^)) which refers for 
fixed index s to the origin of the sth scatterer and not any longer to the origin of 
the 2-dimensional plane. 

Let us explicitly extract the product of the determinants of the subsystems 
from the determinant of the total system (36. 16): 



detS - 



Det (M- /DC) 

Det(M) 
Det (M - /DC) nf=idetS" 
DetM 



^ detS^ 



]~[detS^ 



Det(M-/DC)/nf=idetS 
DetM 



(36.17) 



The final step in the reformulation of the determinant of the S-matrix of the A'^- 
scatterer problem follows from the unitarity of the S-matrix. The unitarity of 
S^(^*) implies for the determinant 



det(S(r)''') - l/detS(yt), 



(36.18) 



where this manipulation is allowed because the T-matrix is trace-class. The uni- 
tarity condition should apply for the S-matrix of the total system, S, as for the 
each of the single subsystems, S% s = I, - ■ ■ ,N. In terms of the result of (36.17), 
this implies 



Det{M{k)-iD{k)C{k)) 
nlidetS-^ 



= Det(M(r)''') 
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since all determinants in (36.17) exist separately and since the determinants det S * 
respect unitarity by themselves. Thus, we finally have 



where all determinants exist separately. 

In summary: We assumed a scattering system of a finite number of non- 
overlapping scatterers which can be of different shape and size, but are all of 
finite extent. We assumed the trace-class character of the T-matrix belonging to 
the total system and of the single-traversal matrix A and finally unitarity of the 
S-matrices of the complete and all subsystems. 

What can one say about the point-particle scattering from a finite number of 
scatterers of ai^bitrary shape and size? As long as each of A'^ < oo single scatterers 
has a finite spatial extent, i.e., can be covered by a finite disk, the total system 
has a finite spatial extent as well. Therefore, it too can be put insided a circular 
domain of finite radius b, e.g., inside a single disk. If the impact parameter of the 
point particle measured with respect to the origin of this disk is larger than the disk 
size (actually larger than (e/2) x b), then the T matrix elements of the A/^-scatterer 
problem become very small. If the wave number k is kept fixed, the modulus of 
the diagonal matrix elements, |r„„„| with the angular momentum m > {e/2)kb, is 
bounded by the corresponding quantity of the covering disk. 

36.2.2 Quantum cycle expansions 

In formula (36. 19) the genuine multi-scattering terms are separated from the single- 
scattering ones. We focus on the multi-scattering terms, i.e., on the ratio of the 
determinants of the multi-scattering matrix M = 1 - A in (36.19), since they are 
the origin of the periodic orbit sums in the semiclassical reduction. The reso- 
nances of the multi-scattering system are given by the zeros of DetM(^) in the 
lower complex wave number plane. 

In order to set up the problem for the semiclassical reduction, we express the 
determinant of the multi-scattering matrix in terms of the traces of the powers 
of the matrix A, by means of the cumulant expansion (36.14). Because of the 
finite number A'^ > 2 of scatterers tr (A" ) receives contributions corresponding to 
all periodic itineraries siS2S3 - ■ ■ s„-is„ of total symbol length n with an alphabet 
Si e {1,2,..., N}. of A^^ symbols. 




(36.19) 



(36.20) 




=-oo =-°° 
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Remember our notation that the trace tr (• • •) refers only to the |/) space. By con- 
struction A describes only scatterer-to-scatterer transitions, so the symbolic dy- 
namics has to respect the no-self-reflection pruning rule: for admissible itineraries 
the successive symbols have to be different. This rule is implemented by the factor 
\-6''' in (36.13). 

The trace tr A" is the sum of all itineraries of length n, 
trA" - trA"'"2A*2^3...A-'„-i^„A*"*i . (36.21) 

{SiS2-S„] 



We will show for the A^-disk problem that these periodic itineraries correspond 
in the semiclassical limit, kog. » 1, to geometrical periodic orbits with the same 
symbolic dynamics. 

For periodic orbits with creeping sections the symbolic alphabet has to be 
extended, see sect. 36.3.1. Furthermore, depending on the geometry, there might 
be nontrivial pruning rules based on the so called ghost orbits, see sect. 36.4.1. 

36.2.3 Symmetry reductions 

The determinants over the multi-scattering matrices run over the multiple index L 
of the multiple index space. This is the proper form for the symmetry reduction 
(in the multiple index space), e.g., if the scatterer configuration is characterized 
by a discrete symmetry group G, we have 

DetM- ]~[(detMD„(^))'''' , 

a 

where the index a runs over all conjugate classes of the symmetry group G and 
Dq. is the ffth representation of dimension da- The symmetry reduction on the 
exact quantum mechanical level is the same as for the classical evolution oper- 
ators spectral determinant factorization (21.17) of sect. 21.4.2. 



36.3 Semiclassical 1-disk scattering 

We start by focusing on the single-scatterer problem. In order to be concrete, we 
will consider the semiclassical reduction of the scattering of a single disk in plane. 

Instead of calculating the semiclassical approximation to the determinant of 
the one-disk system scattering matrix (36.9), we do so for 

A{k) = -^4-lndetSHM = -!-4tr ( InS'(M) (36.22) 
2m dk 2m dk ^ ' 
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the so called time delay. 



_ a -s^ 
2ni ^ 



, H^i^^ka) H^^\ka) 



(36.23) 



Here the prime denotes the derivative with respect to the argument of the Hankel 
functions. Let us introduce the abbreviation 

= ~Q:^ m • (36.24) 

We apply the Watson contour method to (36.23) 

m = y X,n - t: (f) dy -—-Xv ■ (36.25) 
2m ^ 2m 2i Jr sm(y;r) 

m=—oa ^ ^ 

Here the contour C encircles in a counter-clockwise manner a small semiinfinite 
strip D which completely covers the real v-axis but which only has a small finite 
extent into the positive and negative imaginary y direction. The contour C is then 
split up in the path above and below the real y-axis such that 

m = -T-.i dv——yy- dy-— . 

4m tJ-co+,e sm(y7r) J-oo-ie sm(y7r) J 

Then, we perform the substitution v — » -y in the second integral so as to get 

d(^) = t\ \ dy X v + dv . y-y\ 

4n [J_oo+,e sm(y7r) sm(y7r) J 

where we used the fact that x-v = Xv The contour in the last integral can be de- 
formed to pass over the real y-axis since its integrand has no Watson denominator. 

We will now approximate the last expression semiclassically, i.e., under the 
assumption ka :s> 1. As the two contributions in the last line of (36.26) differ by 
the presence or absence of the Watson denominator, they will have to be handled 
semiclassically in different ways: the first will be closed in the upper complex 
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plane and evaluated at the poles ofxv, the second integral will be evaluated on the 
real y-axis under the Debye approximation for Hankel functions. 

We will now work out the first term. The poles of Xv in the upper complex 
plane are given by the zeros of H^y\ka) which will be denoted by V[{ka) and by 

(2) 

the zeros of Hy (ka) which we will denote by -V[{ka), ^ = 1, 2, 3, • • •. In the Airy 
approximation to the Hankel functions they are given by 

V[(ka) - ka + iac{ka) , (36.27) 
-Viika) = -ka + i{at{k*a)y = - {vc{k*a))* , (36.28) 



with 



ia^ika) 



ka 



1/3 



+ e 3 — 



1 



ka 



1/3 ,2 



1 



180 lOka 



30 i 



ka 3150 



29^^ 281^; 



,4 ^ 



62 180 • 63 



(36.29) 



Here qf labels the zeros of the Airy integral 

A(q) = I dr cos(^T - r^) - 3"^/^7rAi(-3"i/^g) , 
Jo 

with Ai(z) being the standard Airy function; approximately, q[ x 6^^^[37t{£ - 
1/4)1^/^/2. In order to keep the notation simple, we will abbreviate Vf = V({ka) 
and Vf = V({ka). Thus the first term of (36.26) becomes finally 



In the second term of (36.26) we will insert the Debye approximations for the 
Hankel functions: 



hI^^^\x) ~ ^ \ — ^ exp (+/ Vx^ - v2 + iv arccos - + /-) for \x\ > v 

(36.30) 



Hi^^^\x) ~ +i J — ^ exp (- Vv^ - + vArcCosh-) for |x| < v . 

V ;r - r2 V xJ 



Note that for v > ka the contributions in Xv cancel. Thus the second integral of 
(36.26) becomes 



2m J_oo 2m J_j,„ a dk \ 



^2 - v2 _ y arccos — + 
kaJ 



1 r^" I 

= dyV^2^2r^ + ... = k + --- , (36.31) 

kn J.ka 2 
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where • • • takes care of the polynomial corrections in the Debye approximation 
and the boundary correction terms in the v integration. 

In summary, the semiclassical approximation to d(^) reads 

d(k) = 2a > — + — \-—k+---. 

Using the definition of the time delay (36.22), we get the following expression for 
detS'CM: 

In det S ' (ka) - Urn In det S ^ (koa) (36.32) 

\ 



. ■ C\r ak 
^ 2ma dk + 2 > 

Jo 2 



e=i 

■k 



+ 



+ 



I — gi2nvc(ka) y _ ^-i2nvc(ka) j 

-lmN{k)+l 2] d^ — [- In + In (l-g-'^'^^^^'-''^)) + • • • 

where in the last expression it has been used that semiclassically ■^V({ka) ~ 
^ve{ka) ~ a and that the Weyl term for a single disk of radius a goes like 
N{k) = na^k^/{4n) + ■ ■ ■ (the next terms come from the boundary terms in the 
v-integration in (36.31)). Note that for the lower limit, ~^ 0> we have two 
simplifications: First, 



1 — Tl\^^ Qcqci) . 
lim 5„,„,,(^0«) ^ lim jr. 5mm' ^ Ix^mm' Vm, m 



— > lim detS'(/:ofl:) - 1 . 
Secondly, for 0, the two terms in the curly bracket of (36.32) cancel. 



36.3.1 1-disk spectrum interpreted; pure creeping 

To summarize: the semiclassical approximation to the determinant S^(^a) is given 
by 

detSkM ~ ; L , (36.33) 



with 



yf(A:a) - ka + iadka) = ka + e'^"^^^{ka/6y^^qe + 
nika) = ka-i(ae{k*a)y = ka + e-'^'^ika/Gy'^qe + 
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Figure 36.2: Right- and left-handed diffractive 
creeping paths of increasing mode number C for 
a single disk. 





and N{ka) = {na^k^)IAn + ■ ■ ■ the leading term in the Weyl approximation for 
the staircase function of the wavenumber eigenvalues in the disk interior. From 
the point of view of the scattering particle, the interior domains of the disks are 
excluded relatively to the free evolution without scattering obstacles. Therefore 
the negative sign in front of the Weyl term. For the same reason, the subleading 
boundary term has here a Neumann structure, although the disks have Dirichlet 
boundary conditions. 

Let us abbreviate the r.h.s. of (36.33) for a disk s as 



detS'(A:a,) ~ (e 



y^s)* 



(36.34) 



where Z^^(kas) and Z^(kas) are the diffractional zeta functions (here and in the fol- 
lowing we will label semiclassical zeta functions with diffractive corrections by a 
tilde) for creeping orbits around the sth disk in the left-handed sense and the right- 
handed sense, respectively (see figure 36.2). The two orientations of the creeping 
orbits are the reason for the exponents 2 in (36.33). Equation (36.33) describes 
the semiclassical approximation to the incoherent part {= the curly bracket on the 
r.h.s.) of the exact expression (36.19) for the case that the scatterers are disks. 

In the following we will discuss the semiclassical resonances in the 1-disk 
scattering problem with Dirichlet boundary conditions, i.e. the so-called shape 
resonances. The quantum mechanical resonances are the poles of the S -matrix in 
the complex ^-plane. As the 1-disk scattering problem is separable, the 5 -matrix 
is already diagonalized in the angular momentum eigenbasis and takes the sim- 
ple form (36.9). The exact quantummechanical poles of the scattering matrix are 
therefore given by the zeros k^^^^ of the Hankel functions H]^\ka) in the lower 
complex k plane which can be labeled by two indices, m and n, where m denotes 
the angular quantum number of the Hankel function and n is a radial quantum 
number. As the Hankel functions have to vanish at specific k values, one cannot 
use the usual Debye approximation as semiclassical approximation for the Hankel 
function, since this approximation only works in case the Hankel function is dom- 
inated by only one saddle. However, for the vanishing of the Hankel function, one 
has to have the interplay of two saddles, thus an Airy approximation is needed as 
in the case of the creeping poles discussed above. The Airy approximation of the 
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Hankel function Hy \ka) of complex-valued index v reads 
n \kaj 

with 

^(^) = e-'M^ {v-ka) + 0[{kar^). 

Hence the zeros of the Hankel function in the complex y plane follow from 
the zeros qc of the Airy integral A{q) (see (36.3). Thus if we set ve = m (with m 
integer), we have the following semiclassical condition on k'''^^ 

m ~ k^'^^'a + iac{k^^''a) 

= 'i—j ik^a] Wo 

[k^^^aj 3150 [ 62 180-63 J 

with / ^ 1,2,3,---. 
For a given index I this is equivalent to 

~ 1 - g(ik"'''-at)2jTa 

the de-Broglie condition on the wave function that encircles the disk. Thus the 
semiclassical resonances of the 1-disk problem are given by the zeros of the fol- 
lowing product 

DO 

1=1 

which is of course nothing else than Zi-disk(^)> the semiclassical diffraction zeta 
function of the 1-disk scattering problem, see (36.34). Note that this expression 
includes just the pure creeping contribution and no genuine geometrical parts. 
Because of 

H^^lika) = {-iy"Hl^\ka), 

the zeros are doubly degenerate if m t 0, corresponding to right- and left handed 
creeping turns. The case m = is unphysical, since all zeros of the Hankel func- 
tion tf^\kd) have negative real value. 



-— f>-^] 

lOk'^^'a \ 30 j 

' ' ' » 

(36.35) 
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Figure 36.3: The shape resonances of the 1-disk 
system in the complex k plane in units of the disk 
radius a. The boxes label the exact quantum me- 
chanical resonances (given by the zeros of H',^\ka) 
for m = 0, 1,2), the crosses label the diffractional 
semiclassical resonances (given by the zeros of 
the creeping formula in the Airy approximation 
(36.35) up to the order Odkay^)). 



Figure 36.4: Same as in figure 36.3. However, 
the subleading terms in the Airy approximation 
(36.35) are taken into account up to the order 
0([kay''^) (upper panel) and up to order (9([fca]"') 
(lower panel). 



From figure 36.3 one notes that the creeping terms in the Airy order 0{[ka]^^^), 
which are used in the Keller construction, systematically underestimate the magni- 
tude of the imaginary parts of the exact data. However, the creeping data become 
better for increasing Rek and decreasing \lmk\, as they should as semiclassical 
approximations. 

In the upper panel of figure 36.4 one sees the change, when the next order 
in the Airy approximation (36.35) is taken into account. The approximation is 
nearly perfect, especially for the leading row of resonances. The second Airy 
approximation using (36.35) up to order 0{[kaY^) is perfect up to the drawing 
scale of figure 36.4 (lower panel). 
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Figure 36.5: A 4-disk problem with three specular 
reflections, one ghost tunneling, and distinct creeping 
segments from which all associated creeping paths can 
be constructed. 




36.4 From quantum cycle to semiclassical cycle 



The procedure for the semiclassical approximation of a general periodic itinerary 
(36.20) of length n is somewhat laborious, and we will only sketch the procedure 
here. It follows, in fact, rather closely the methods developed for the semiclassical 
reduction of the determinant of the 1-disk system. 

The quantum cycle 



trA 



SlS2 



A 



■«»!.« I 



A Si"! 
■'I -'2 



■A 



stiU has the structure of a "multi-trace" with respect to angular momentum. 

Each of the sums YjT =-oa ~ '^^e 1-disk case - is replaced by a Watson 
contour resummation in terms of complex angular momentum y^ .. Then the paths 
below the real y.v,-axes are transformed to paths above these axes, and the integrals 
split into expressions with and without an explicit Watson sin(y v,7r) denominator. 

1. In the sin(yv,7r) -independent integrals we replace all Hankel and Bessel 
functions by Debye approximations. Then we evaluate the expression in 
the saddle point approximation: either left or right specular reflection at 
disk Si or ghost tunneling through disk st result. 

2. For the sin(yi.7r) -dependent integrals, we close the contour in the upper y^,. 
plane and evaluate the integral at the residua //y]\^a.5,)=0. Then we use 

the Airy approximation for Jy^Xkag.) and //yj\^a.v,): left and right creeping 
paths around disk Si result. 



In the above we have assumed that no grazing geometrical paths appear. If 
they do show up, the analysis has to be extended to the case of coinciding saddles 
between the geometrical paths with n/2 angle reflection from the disk surface and 
paths with direct ghost tunneling through the disk. 

There ai^e three possibilities of "semiclassical" contact of the point particle 
with the disk sf. 

1. either geometrical which in turn splits into three alternatives 
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(a) specular reflection to the right, 

(b) specular reflection to the left, 

(c) or 'ghost tunneling' where the latter induce the nontrivial pruning 
rules (as discussed above) 

2. or right-handed creeping turns 

3. or left-handed creeping turns, 

see figure 36.5. The specular reflection to the right is linked to left-handed creep- 
ing paths with at least one knot. The specular reflection to the left matches a 
right-handed creeping paths with at least one knot, whereas the shortest left- and 
right-handed creeping paths in the ghost tunneling case are topologically trivial. 
In fact, the topology of the creeping paths encodes the choice between the three 
alternatives for the geometrical contact with the disk. This is the case for the 
simple reason that creeping sections have to be positive definite in length: the 
creeping ampUtude has to decrease during the creeping process, as tangential rays 
are constantly emitted. In mathematical terms, it means that the creeping angle 
has to be positive. Thus, the positivity of the two creeping angles for the shortest 
left and right turn uniquely specifies the topology of the creeping sections which 
in turn specifies which of the three alternatives, either specular reflection to the 
right or to the left or straight "ghost" tunneling through disk j, is realized for the 
semiclassical geometrical path. Hence, the existence of a unique saddle point is 
guaranteed. 

In order to be concrete, we will restrict ourselves in the following to the scat- 
tering from N < oo non-overlapping disks fixed in the 2-dimensional plane. The 
semiclassical approximation of the periodic itinerary 

tr A'^'^^A*^"*' • • . A^""''*" A*"*' 

becomes a standard periodic orbit labeled by the symbol sequence s\S2 - ■ ■ s„ . De- 
pending on the geometry, the individual legs Si-^i — > Si — > si+i result either from a 
standard specular reflection at disk Sj or from a ghost path passing straight through 
disk Si. If furthermore creeping contributions are taken into account, the symbolic 
dynamics has to be generalized from single-letter symbols {5,} to triple-letter sym- 
bols {si,o-i X with ii > 1 integer valued and cr,- = 0, ±1 ^ By definition, the 
value cTj = represents the non-creeping case, such that {si, x {j} = {si, Oj = {si] 
reduces to the old single-letter symbol. The magnitude of a nonzero corre- 
sponds to creeping sections of mode number whereas the sign cr, = ±1 signals 
whether the creeping path turns around the disk 5, in the positive or negative sense. 
Additional full creeping turns around a disk Sj can be summed up as a geometrical 
series; therefore they do not lead to the introduction of a further symbol. 
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Figure 36.6: (a) The ghost itinerary (1, 2, 3, 4). (b) 
The parent itinerary (1, 3, 4). 




36.4.1 Ghost contributions 

An itinerary with a semiclassical ghost section at, say, disk si can be shown to 
have the same weight as the corresponding itinerary without the Sj th symbol. 
Thus, semiclassically, they cancel each other in the tr ln(l - A) expansion, where 
they are multiplied by the permutation factor n/r with the integer r counting the 
repeats. For example, let (1,2,3,4) be a non-repeated periodic itinerary with a 
ghost section at disk 2 steming from the 4th-order trace trA^. By convention, 
an underlined disk index signals a ghost passage (as in figure 36.6a), with cor- 
responding semiclassical ghost traversal matrices also underlined, a''— A— 
Then its semiclassical, geometrical contribution to tr ln(l - A) cancels exactly 
against the one of its "parent" itinerary (1,3,4) (see figure 36.6b) resulting from 
the 3rd-order trace: 

_1 (4A1'^A^-3a3'4a4'1) - 1 (3 a1'3^3,4a4,iJ 
= (+1-1)A''3a3-4a4'1 =0. 

The prefactors -1/3 and -1/4 are due to the expansion of the logarithm, the fac- 
tors 3 and 4 inside the brackets result from the cyclic permutation of the periodic 
itineraries, and the cancellation stems from the rule 

. . . A'-i±iA^''+2 ... = ... (-A'''+2) .... (36.36) 

The reader might study more compUcated examples and convince herself that the 
rule (36.36).is sufficient to cancel any primary or repeated periodic orbit with one 
or more ghost sections completely out of the expansion of tr ln(l - A) and there- 
fore also out of the cumulant expansion in the semiclassical limit: Any periodic 
orbit of length m with n{< m) ghost sections is cancelled by the sum of all 'par- 
ent' periodic orbits of length m - i (with \ < i <n and / ghost sections removed) 
weighted by their cyclic permutation factor and by the prefactor resulting from 
the trace-log expansion. This is the way in which the nontrivial pruning for the 
A^-disk billiards can be derived from the exact quantum mechanical expressions 
in the semiclassical limit. Note that there must exist at least one index / in any 
given periodic itinerary which corresponds to a non-ghost section, since other- 
wise the itinerary in the semiclassical limit could only be straight and therefore 
nonperiodic. Furthermore, the series in the ghost cancelation has to stop at the 
2nd-order trace, tr A^, as tr A itself vanishes identically in the full domain which 
is considered here. 

'Actually, these are double-letter symbols as cr, and Z, are only counted as a product. 
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36.5 Heisenberg uncertainty 

Where is the boundary ka w 2'^~^L/a coming from? 

This boundary follows from a combination of the uncertainty principle with 
ray optics and the non-vanishing value for the topological entropy of the 3-disk 
repeller. When the wave number k is fixed, quantum mechanics can only resolve 
the classical repelling set up to the critical topological order n.The quantum wave 
packet which explores the repelling set has to disentangle 2" different sections 
of size d ~ a/2" on the "visible" part of the disk surface (which is of order a) 
between any two successive disk collisions. Successive collisions are separated 
spatially by the mean flight length L, and the flux spreads with a factor L/a. In 
other words, the uncertainty principle bounds the maximal sensible truncation in 
the cycle expansion order by the highest quantum resolution attainable for a given 
wavenumber k. 



Commentary 

Remark 36.1 Sources. This chapter is based in its entirety on ref. [K. f]; the reader 
is referred to the full exposition for the proofs and discussion of details omitted here, 
sect. 36.3 is based on appendix E of ref. [K. f]. We follow Franz [35.19] in applying the 
Watson contour method [35.20] to (36.23). The Airy and Debye approximations to the 
Hankel functions are given in ref. [35.21], the Airy expansion of the 1-disk zeros can be 
found in ref [35.22].For details see ref s. [35.19, 35.22, 35.23, K.l]. That the interior do- 
mains of the disks are excluded relatively to the free evolution without scattering obstacles 
was noted inrefs. [35.24, 35.15]. 

The procedure for the semiclassical approximation of a general periodic itinerary 
(36.20) of length « can be found in ref. [K.l] for the case of the A^-disk systems. The 
reader interested in the details of the semiclassical reduction is advised to consult this 
reference. 

The ghost orbits were introduced in refs. [35.12, 35.24]. 

Remark 36.2 Krein-Friedel-Lloyd formula. In the literature (see, e.g., refs. [35.14, 
35.15] based on ref. [35.11] or ref. [35.1]) the transition from the quantum mechan- 
ics to the semiclassics of scattering problems has been performed via the semiclassical 
limit of the left hand sides of the Krein-Friedel-Lloyd sum for the (integrated) spectral 
density [K.5, K.6, 35.8, 35.9]. See also ref. [35.13] for a modern discussion of the ICrein- 
Friedel-Lloyd formula and refs. [35.1, 35.17] for the connection of (35.17) to the Wigner 
time delay. 

The order of the two limits in (35.18) and (35.17) is essential, see e.g. Balian and 
Bloch [35.1 1] who stress that smoothed level densities should be inserted into the Friedel 
sums. 

The necessity of the +;e in the semiclassical calculation can be understood by purely 
phenomenological considerations; Without the ie term there is no reason why one should 
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be able to neglect spurious periodic orbits which solely are there because of the introduc- 
tion of the confining boundary. The subtraction of the second (empty) reference system 
helps just in the removal of those spurious periodic orbits which never encounter the scat- 
tering region. The ones that do would still survive the first limit ^ oo, if they were not 
damped out by the +;e term. exercise 35.1 



Remark 36.3 T, C, D' and A"' matrices are trace-class In refs. [K.l] it has 
explicitly been shown that the T-matrix as well as the C, D' and A**'-matrices of the 
scattering problem from N < co non-overlapping finite disks are all trace-class. The 
corresponding properties for the single-disk systems is particulary easy to prove. 
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